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ABSTRACT 



We exploit our formula for the gravitational potential of finite size, power-law disks to derive a general expression 
linking the mass of the black hole in active galactic nuclei (AGN), the mass of the surrounding disk, its surface density 
profile (through the power index s), and the differential rotation law. We find that the global rotation curve v(R) of the 
disk in centrifugal balance does not obey a power law of the cylindrical radius R (except in the confusing case s = — 2 
that mimics a Keplerian motion), and discuss the local velocity index. This formula can help to understand how, from 
position-velocity diagrams, mass is shared between the disk and the black hole. To this purpose, we have checked the 
idea by generating a sample of synthetic data with different levels of Gaussian noise, added in radius. It turns out 
that, when observations are spread over a large radial domain and exhibit low dispersion (standard deviation a < 10% 
typically), the disk properties (mass and s-parameter) and black hole mass can be deduced from a non linear fit of 
kinematic data plotted on a (R, 7?u 2 )-diagram. For a > 10%, masses are estimated fairly well from a linear regression 
(corresponding to the zeroth-order treatment of the formula), but the power index s is no longer accessible. We have 
applied the model to 7 AGN disks whose rotation has already been probed through water maser emission. For NGC3393 
and UGC3789, the masses seem well constrained through the linear approach. For IC1481, the power-law exponent s 
can even be deduced. Because the model is scale-free, it applies to any kind of star/disk system. Extension to disks 
around young stars showing deviation from Keplerian motion is thus straightforward. 

Key words. Accretion, accretion disks — Gravitation — Methods: analytical — Galaxies: active 



1. Introduction 



up t o instability (e.g. 



The mass of astrophysical objects — except maybe for stars 
— is generally difficult to determine with precision, mostly 
because of inappropriate tracers, relatively low spatial res- 
olution, and a certain misunderstanding of the internal 

structure and physical processes involved. This is the case |lQ87l : iPapaloizou fc Lin 
for giant disks orbiting supermassive black holes in active 
galactic nuclei (AGN). For some nearby objects, the cold 
gas rotating in the outermost regions (the subparsec scale 
typically from the center) is detected at radio wavelengths 



The mass of disks can also be probed via the conse- 
quences of their gravity. All the material contained within 
a disk exerts gravitational forces on itself — the so-called 
"self- gravity" — which influences or even strongly gov- 
erns (like in galaxies) i nternal orbita l motions, sometimes 



Mestell Il963t iBinnev k, Tremaind 
119951) . Even in the presence of 



through water vap or emission (e.g. iMivoshi et all 1 19951: 
iBraatz et all 12009). The inner regions, not accessible yet 



to current instruments, could host the bulk of the mass if 
the total surface density in the disk va ries roughly with the 
cylindrical radius R as R~ 2 or faster (|Shakura fc Sunvaevl 
ll973tlCollin-Souffrin fc Dumontlll99"ol ). Estimating the disk 
mass is a complex task. It necessitates a global disk model 
capable of describing the dynamics of gas, its thermody- 
namics, its chemical complexity, as well as its interaction 
with radiation (lines and continuum). The disk mass is 
an important quantity in understanding the AGN phe- 
nomenon. Along with the accretion rate, turbulent viscos- 
ity, and black hole mass, it helps to put constraints on the 
activity of the AGN in terms of stability, lifeti me, lumi- 



nosity, and matter supplie d from the host galaxy ()Combesl 
120011: iGollin fc Zahdl2008lh 
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a massive central object, self- gravity may cause a slight 
deviation in Kepler's law, which is interesting to ana- 
lyze and to quantify. Obviously, non-Keplerian rotation 
can have other o rigins like pressur e effec ts such as in 
slim/t hick disks (lAbramowicz etaLl EUl) or magnetic 
fields ([Hevvaerts k Priest! 19891 ). Here, we focus on self- 
gravity, whi ch is expected to play a role in geometrically 
thin disks (|Shore k White! 119821 IShlosman fc Begelmanl 
119871) . 

There are many articles that aims to establish the re- 
lation between the dynamics (through gravitational po- 
tentials) and mass density dis tribution, especially in th e 
context of galactic dynamics (Bi nnev fc Tremaind Il987). 
Existing potential/density pairs do not however seem fully 
appropriate to gaseous disks surrounding a central ob- 
ject, probably because star/disk systems where the gas ex- 
hibit non-Keplerian mot ions are still marginal. As models 
and theories suggest (e.g-FShakura k Sunvaevlll973l ; |Pringld 
119811: iGollin-Souffrin fc Dumontlll990l lHurelll998lV gaseous 
disks in AGN are large and are expecte d to exhibit aself- 
similar behavior over some radial range. IHure eTaTl ((20081) 
has determined an accurate formula for the gravitational 
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potential in the midplane of a flat power-law disk with fi- 
nite size and mass. It is valid for a wide range of the power 
index for surface density. In this article, we use this result 
to derive an algebraic relation between the orbital veloc- 
ity of the gas, the disk parameters (surface density profile, 
mass, size), and the mass of the central object, assuming a 
pressure-less disk at centrifugal equilibrium, as commonly 
done. This relation furnishes a simple method for deter- 
mining how the mass is shared between the disk and the 
central object. As expected, the "modified" rotation law is 
not, in this model, a po wer law of the radius a s often con- 
sidered in this context (|Herrnstein et all |2005|) . Although 
this study is valid for any kind of astrophysical star/disk 
system (like in circumstellar environments), we focus on 
AGN disks whose kinematics have been observed in VLBI 
through water-maser emission. 

This article is organized as follows. In Section [21 we re- 
call the model of a pressure-less disk at centrifugal equilib- 
rium surrounding a central black hole. We introduce the 
formula for the gravitatio nal potential in the disk mid- 
plane bv lHure et al.l ((2008), and derive the general expres- 
sion for the velocity of the orbiting g function of 
the disk mass, black hole mass, and surface density profile 
through the power-law index. In Section [3j we show how 
this expression (or its zero-order version) can be used to 
estimate how the mass is shared between a central object 
and its surrounding disk by analyzing observational data 
in a "position-dynamical mass" diagram (instead of the 
classical position- velocity diagram). We first applied the 
method t o IC1481, thereby refini ng the disk parameters re- 
ported in iMamvoda et al.l (|2009[) . We discuss uncertainties 
in Sect. HJ and show how dispersion naturally goes against 
the method. Section [5] is devoted to applying of the method 
to a sample of a few well known AGN hosting a masing 
outer disk. We conclude in Section H3 



2. The basic model 

2.1. A pressure-less disk at centrifugal equilibrium 

We consider a gaseous disk with inner edge ai n and outer 
edge a out 3> &i n , orbiting a central black hole with mass 
-AIbh- This disk is assumed to be axially symmetrical, flat 
(i.e. no vertical thickness), pressure-less, and steady. At cen- 
trifugal equilibrium, the rotation velocity v of material at 
cylindrical distance R in the midplane of the disk, in the 
reference frame of the black hole, is given by the standard 
relation: 



GM] 



BH 



R 



R 



dipd 
dR 



(1) 



where ipd is the gravitational potential of the disk. This lat- 
ter function critically depends on the surface density pro- 
file S(i?) through the Poisson integral. It is generally not 
easily accessed by analytical means, even in the actual one- 
dimensional case. 

There is a broad literature devoted to determin- 
ing potential- density pairs {il)&,Yi) for axially symmetric 
syste ms (e.g. iBinnev fc TremaineT Il987t lEvans fc Collettj 
[1991 . Here, we consider the class of flat, power-law dis- 
tributions where the surface density varies according to 
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Fig. 1. P2{s) in the range of validity of Eq. [3J 



where R — a ou tzu, A = 0i n /a O ut is the axis ratio, E out 
the surface density at the outer edge, and s is a constant. 
Such a profile seems well-suited for large, gaseous disks in 
AGN, at least in the framework of ge ometrically thin disk 
models that predict s ~ — 1 typically (IShakura fc Sunvaev 
19731 lPrinridll98lt ICollin-Souffrin fc Dumontl 119901: IHure 



1998). Potential-density pairs for flat power-law disks, in- 
cluding Mestel ' s solu tion (s = —1), are summarized in 
lEvans fc Read! (|1998f) . These correspond to infinite disks 
(i.e. cii n = a nd o ou t — > 00) whose mass is infinite as 
soon as s > —2. IConw av (2000) has produced formal solu- 
tions corresponding to finite disks and no inner edge (i.e. 
«in = 0), but for even positive indexes (i.e. s = 0,2,...). 
Unfortunately, acco unting for edges i ncreases the mathe- 
matical difficulties. IHure et ail ()2008l ) have recently pro- 
duced a reliable approximation for -0d associated with Eq. 
12 namely (see their Eq. 53) 



1 w 



2+s 



A 2+s 1 - w 1+s 



wherf] Vw = 27rGS 0Ut a ut and 



B = 6C — 1 « 0.431, 



(3) 



(4) 



where C is the Catalan's constant. This approximation is 
accurate within a few percent, provided the disk is large 
enough (A -C 1) and —3 < s < 0. These condi- 
tions are probably me t in most astrophysical d i sks, es- 
pecially in AGN disks (IShakura, fc Sunvaevlll973t iPringld 
ll98lUCollin-Souffrin fc Dumondll990UHure1ll998lV and oth- 
ers (e.g. lDubrulldll992l ). Actually, we have A ~ 3i?s / a ou t ~ 
10~ 5 for a 10 8 M Q AGN black hole accreting a parsec size 
disk (Rg being the Schwarzschild radius). Besides, Eq. [3] 
does not suffer from the edge singularities expected when 
considering sharp edges. 



elsewhere, 



[A,l], 



(2) 



1 This constant ip ou t has dimension of a potential, but it is not 
the value at the outer edge. 
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2.2. The dynamical mass 

As Eq. [JJ shows, a "good" variable to m easure masses in 
this model is (e.g. lYamauchi et al.ll2004l) : 



Rv - ( \ 



(5) 



and the disk makes its own contribution through -0d- In 
the case of a spherical distribution (where R —> r), this 
dynamical mass would represent t he enclosed mass at a 
given radius (|Herrnstein et al.]|2005l ). Here, things are quite 
different since matter is gathered in a plane. There is no 
obvious use of the Gauss theorem, although the monopole 
approximation gives the right order of magnitude. Inserting 
in Eq. [5] the velocity v given by Eq. [TJ and using Eq. [3] for 
the disk potential ipd, we finally get 



H{m) = M BH - M< 



A 2+s -m s+2 P 2 {s) 
1 - A 2 + s ' 



(G) 



where Md is the disk mass (at a ou t), and 

P 2 {s) = l-B(l + s)(2 + s) (7) 

is a second-order polynomial in s. It is displayed in Fig. [JJ 
In the range of interest, P 2 (s) is always positive and has a 
limited range of variation, since 0.14 < P 2 (s) < 1.11. We 
also have P 2 (— 2) = P2(— 1) = 1 and a mean value of ~ 0.8. 

From Eq. [UJ we expect two extreme behaviors of the 
function /j.(zj) depending on the disk properties. If s < —2 
(the case of "centrally peaked" distributions), we have 

(8) 



Mi 



BH 



est. 



In the absence of radial gradient of /x, it is not possible 
to separate the disk and the central black hole, regardless 
of the disk mass. A Keplerian rotation curve result from a 
massive disk (without central object) or to a massive central 
object (with a light disk, as often considered). If s > — 2 
(the bulk of the disk mass stands in the outermost regions) , 
we have 



~ M BH + M d m s+2 P 2 (s), 



(9) 



which is essentially an increasing function of the radius 
(positive gradient). Disk rotation is therefore always super- 
Keplerian. The larger the disk mass, the larger the devia- 
tion from Kepler's law. 

We conclude that, if the disk mass is significant with 
respect to the central mass, v 2 (R) is the sum of two power 
laws of the radius, and this sum is not a power law (see 
below). In other words, fitting the global rotation curve of 
a system containing a massive disk and a black hole with a 
single power law (e.g. v oc i? 7 ), cannot give any quantitative 
information about the mass distribution in the framework 
of Newtonian gravity. This approach is often considered 
when a massive disk is suspe c ted (see e.g. Herrnstein et al.l 
l2005t iKondratko et~all[200l iMcCallum et all 12009ft . S~ 
however, appendix [X] for a short discussion of the velocity 
index 7. 

For s = —1, Eq. Oreads 



H{w) = M BH + Mi 



w - A 



(10) 



which is to be compared to the case of a Mestel disk (jMestell 
[1961 

AtMestcl(n7) = M' A {w), (11) 



where M d is the cumulative disk mass at the actual ra- 
dius (linea r with the radius w). As already pointed out 
elsewhere (|Binnev fc Tremaind Il987t ) . this equation "hap- 
pens to give the same answer" as what is deduced from the 
Gauss theorem for a spherically symmetric distribution. In 
contrast, Eq. [TUJ corresponds to a Mestel disk truncated on 
both sides. It includes edge effects and total disk mass, and 
explicitly contains the central point mass (which is not part 
of Mestel's disk model). In the limit A — > 0, we then recover 
Mestel's solution. 



3. AGN disk/black hole mass determination 
method. The case of IC1481 

3.1. Position-dynamical mass diagram 

We immediately see from above that, if the rotation curve 
of the disk is partly known in the form of N observa- 
tional points {(Ri,Vi)}, then some constraints can be set 
on the disk mass, black hole mass, and surface density 
profile by fitting the data through Eqs. [5] or 

O Obviously, this procedure does not guarantee that the 
triplet (Mbh 7 Md, s) is physically meaningful given the sim- 
plicity of our model and assumptions. Uncertainties in data 
also fragilize the inversion. Thus, there are three different 
possibilities. 

A: Data points {(Ri, Hi)}n show no noticeable variation 
around a constant value, only a certain dispersion. The 
systems thus appears in Keplerian rotation. We deduce 
that either there is a light disk surrounding a massive 
black hole or the disk is rather massive but the gas is 
distributed such that s < — 2. The diagram only gives 
the quantity Mbh+Mj, which can be identified with the 
so-called "binding mass" Mb (or enclosed mass). There 
is no way to separate the black hole and the disk in this 
analysis. 

B: Data points show a significant variation, still with a 
certain dispersion (see Sect. [4]). The gas rotates faster 
than Keplerian. If fi increases faster than zu, then s > 
— 1, otherwise s < —1. In either case, fitting the data 
points through Eq. [9] can yield a triplet (Mbh, Md, s). 

C: Data cannot be fitted by Eq. or inferred parameters 
are non physical. In this case, our model is inappropri- 
ate. Various reasons can be invoked (see Sect. [6]). 

3.2. Zeroth-order: disk mass and central mass 

The zeroth-order treatment of the non linear formula is 
interesting and instructive because it gives the orders of 
magnitude. Actually, if we consider that astrophysical disks 
are characterized by s ~ — 1, we can expand Eq. [S] around 
s = -1. We find (see also Eq. [TU] with A -> 0) 



M BH 
M bh 



M d vj - 
Mdtu. 



Mi(s + 1)zu\d.zu, 



(12) 



We conclude that, if observational data plotted on a 
position-dynamical mass diagram are almost linearly dis- 
tributed in a position-dynamical mass diagram, then the 
slope is the disk mass Mi and the intercept is the black 
hole mass. In the following, we discuss both approaches in 
the context of AGN disks whose rotation, for some of them, 
is known from maser emission. 
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Fig. 2. Position- velocity data of maser spots in the disk of 
IC1481 (top), and corresponding position-dynamical mass 
diagram (bottom). A non linear fit of the data (green 
curve) gives the disk mass, the black hole mass and the 
s-param eter. The linear f i t (blac k line) is also shown. Data 
are from iMamvoda et al.l (1200S 



3.3. Scale-free formula. Scaling to AGN disks 

Equation [5] is totally scale- free if we divide fj, by the black 
hole mass. Actually, if q — M^/Mbh denotes the mass ratio 
in the system, Eq. [9] writes 



li(w) 



M 



BH 



l + qm s+2 P 2 (s), 



(13) 



so that the model applies to any kind of star /disk system. In 
the context of AGN, black holes are supermassive (several 
million solar masses), the central accretion disk typically 
has the parsec size, and rotational velocities are, at such 
distances, hundreds of km/s. With this scaling, the formula 
for the dynamical mass become^] 



Mi 



10 7 M, 



o 



0.232 



Ri 
1 pc 



100 km/s 



(15) 
(16) 



where Vi is, following our model, the disk rotational velocity 
measured relative to the systemic/receding velocity of the 
system, and Ri relative to the rotation axis of the disk. In 



For a circumstellar system, we would have, for instance, 

2 



/J; 



\M, 



1.12 







Ri 



100 AU J V 100 km/s 



(14) 



AGN points (m) M B h M d 



IC1481 
IC1481 



26 
26 



3.67 
3.67 



1.51 
1.30 



3.59 
3.55 



-0.88 
(-1) 



0.89 
0.89 



Table 1. Results of non linear fit through Eq. [17] for 
IC1481 (top), and results for the linear fit (a least-square 
fit) through Eq.[T5]( bottom). La st column is the correlat ion 
coefficient. Maser data are fromlMamvo da et al.l (J2009). 



the following, masses are expressed in units of 10 7 M Q , and 
denoted M. Then, Eq. becomes 

AM = M BH + M d m s+2 P 2 (s), 



(17) 
(18) 



and its linear version is 

jl(w) = M B h + M d ru. 

3.4. An example: the case of IC1481 

The top panel in Fig. [5] shows the rotational velocity of 
maser spots observed in the outer disk of the active nu- 
cleus in galaxy IC1481 versus the distance from the center 
(jMamvoda et al.|[20 09) . The radius of the outermost maser 
spot is set to a out (although the gas disk can extend farther 
awajQ) ■ The bottom panel displays the same data once con- 
verted into a (tn, /i)— diagram, as well as the result of the 
fit by Eq. [17] (non-linear) and by Eq. [18] (linear). For each 
point, the square of the radial deviation between the data 
and the linear fit is shown (see below). The parameters of 
the two fits are gathered in Table [T] These are in good 
agreement, especially because the solution of the non lin- 
ear fit gives an s-parameter close to — 1. In such a case, 
the curvature is difficu lt to detect by eye. As announced in 
IMamvoda et~aTl (|2009lffl . the disk mass is higher than the 
black hole mass by a factor of 2 — 3. 

4. Uncertainties and data dispersion 

Position- velocity data deduced from observations generally 
suffer from dispersion and uncertainties, which can have 
different origins: physical (e.g. variability, non uniform dy- 
namics, geometry and deprojection) and instrumental (i.e. 
lack of resolution). In particular, locating t he position of 
emitt ers precisely is a critical point (e.g. lUscanga et al.l 
l2007f ). In order to check whether the method is "robust" 
enough to infer some reliable information about masses, 
we generated a sample of N points {tUi,fli}, obeying ex- 
actly Eq. [T7]for a given reference triplet (Mbh, Md.) s) ref . 
To these synthetic data, we added uncertainties on the ra- 
dius (without presuming their origin) with three different 
levels of "noise" . We considered a Gaussian noise, with var- 
ious standard deviations a. The perturbed radii roj were 
then all rescaled so that finally tujv = 1. Unfortunately, 
this rescaling procedure introduce a slight bias by tending 
to overestimate both the black hole mass mass and the disk 



3 It can be shown that Eqs. \E\ and [5] still hold if one uses a 
reference value an < a out other than the outer radius a out since 
Md oc R a+2 . Then M<j refers to the cumulative disk mass up to 
this reference radius an- 

4 In IMamvoda et all {2009), an extreme case without black 
hole was considered. 
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Fig. 3. Synthetic data set made with N = 20 points gen- 
erated from Eq. 13 with (Af B H, M d , s) rcf = (1,4,-1.5) 

(top) and with (Af B H, M d , s) rcf = (4,1,-1.5) {bottom). 
Uncertainties corresponding to a Gaussian noise with a 
standard deviation cr = 0.03 have been added to all the 
radii. Data were rescaled in radius such that the outer point 
has a normalized radius zun = I. 



Here, we report a typical experiment obtained for 



a massive disk configuration with (Mbh> Md) rof 
corresponding to mass ratio q — 4, 
a massive black hole configuration with (Mbh> Md) Icl 
(4, 1) corresponding to q = i, 



(1,4) 



with s rcf = —1.5 and N — 20 points in both cases. Initially, 
data were randomly spread over the range [0.5, 1], and we 
considered 3 levels of noise: a = {0.03, 0.05, 0.1}. The syn- 
thetic (w 1 /i)-diagram obtained for a = 0.03, which is pre- 
sumed to mimic dispersed observationnal data, is shown in 
Fig. [3] with the exact rotation law given by Eq.[17] We tried 
to determine the best triplet (M B n,Md,s) in three ways: 



— by the linear approach, 

— by the non linear approach, but forcing s 

— by the non linear approach. 



ref 



The second case is only illustrative, as it can not be used 
in practice since s is not known a priori. For each config- 
uration, the relative distance between the reference triplet 
and the one deduced by fitting the synthetic noisy data is 
measured by the parameter 8, with 



AL 



BH 



Md 



As a result, <5 — ^ when the method works, whereas 
8 > 1 when it fails (i.e. when the output and the in- 



a 


0.03 


0.05 


0.1 


comments 


S 


1.79 


1.67 


1.24 


linear approach 


Mbh 


2.74 


2.62 


2.19 


(s = -1) 


Md 


2.89 


3.11 


3.91 




5 


0.33 


0.61 


1.57 


non-linear 


M BH 


0.68 


0.42 


-0.49 


approach 


M d 


4.43 


4.75 


5.89 


(with s — — 1.5) 


6 


0.32 


0.83 


no solution 


non-linear 


Mbh 


1.31 


0.20 




approach 


M d 


3.88 


4.95 






s 


-1.41 


-1.52 






Massive black hole configuration (q = 




a 


0.03 


0.05 


0.1 


comments 


8 


0.34 


0.46 


1.44 


linear approach 


M BH 


4.28 


4.06 


3.38 


= -i) 


M d 


0.96 


1.31 


2.39 




8 


0.49 


1.03 


2.70 


non-linear 


M BH 


3.59 


3.13 


1.71 


approach 


M d 


1.48 


2.01 


3.64 


(with s = —1.5) 


8 


0.20 


4.52 


no solution 


non-linear 


M BH 


4.07 


-0.45 




approach 


M d 


1.09 


5.46 






s 


-1.24 


-1.83 







Table 2. Results for synthetic data generated from Eq. 
1171 for the massive disk configuration (top), and for the 
massive black hole configuration (bottom). Values in bold 
correspond to the best approach. 



put triplets differ significantlj0) . Practically, the linear ap- 
proach is performed from a least-square procedure. The re- 
sults are gathered in Table [2j As expected, the lower the 
dispersion, the better the method. For the massive disk 
configuration and a — 0.03 (low dispersion), we recover 
the disk mass within a few percent by the non linear fit, 
while the error in the black hole mass is about 30%. By the 
linear approach, we get the disk mass within 30%, while 
there is a factor 2 — 3 for the black hole mass. This is 
not all that surprising since the curvature of the function 
jj.(w) is noticeable in this case, so the error on the inter- 
cept can be large. As the dispersion increases, the linear 
approach becomes more and more reliable for getting the 
disk mass (within a few percent), but the black hole mass is 
still poorly determined (a factor 2 typically). Regarding the 
massive black hole configuration and a = 0.03 (low disper- 
sion), the black hole mass and the disk mass are deduced 
correctly within a few percent typically, especially using 
the non linear method. As the dispersion increases, the un- 
certainty in both two quantities rises, for the same reason 
as mentioned hereabove; however, the black hole mass is 
determined with about 15%. 



J This norm can introduce a bias into the interpretation as, a 
relative discrepancy of 60% in each direction suffices to produce 
5=1, but 8 can be large even if two of the three parameters are 
correct. 
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Table 3. List of AGN hosting a masing disk considered here, references for position-velocity data and associated 
disk/black hole parameters reported by authors model only). Columns 3-5 gives values found by the authors. Masses 
are given in units of 10 7 M Q ; Mb is the binding mass as determined from spherically symmetric models. 



AGN N (jli) Mbh cor. case* st. dev. a comment for non linear fitting 



IC1481 


26 


3.67 


1.30 


3.55 


0.89 


B 


0.11 


agreement 




UGC3789 


39 


1.09 


0.81 


0.62 


0.89 


B 


0.07 


(M B n,M d ,s) p 


a (5.91,-2.13,-0.056) 


NGC3393 


17 


2.85 


0.59 


3.48 


0.95 


B 


0.06 


no solution 




NGC4258 


50 


3.78 


3.63 


0.16 


0.37 


AB 


0.3 


{Mbh, Md, s) s= 


s (3.20,0.60,-1.97), cor. 0.38 


NGC1068 


32 


1.63 


0.72 


1.09 


0.49 


BC 


0.2 


(M B H,M d ,s) p 


s (2.77, -0.89, -3.58) 


NGC4945 


6 


0.09 


0.11 


-0.02 


-0.11 


C 


2 


{Mbh , Md , s ) s= 


s (0.025,0.053,-3.38), cor. 0.14 (s < -2) 


Circinus 


21 


0.12 


0.04 


0.14 


0.37 


c 


0.6 


no solution 





Table 4. Results of linear-fitting through Eq. [TH1 Column 8 gives the standard deviation a. *For the classification, see 
Sect. 1 



Since the errors on the radius propagate to the variable 
/it, we made the same fit in the (i?,i>)-diagram and found 
very similar results. In brief, we see that when the "global 
slope" is steep (massive disk configuration), the disk mass 
is determined with precision. Conversely, when the "global 
slope" is almost zero (massive black hole configuration), 
the black hole mass (value at the intercept) is well deter- 
mined. When dispersion is too large, the correlation be- 
comes too weak for interpretation, and the method fails to 
give the triplet, but it is important to note that the linear 
approach always gives the value of the most massive com- 
ponent within a factor less than about 15%, and the less 
massive by a factor 2 — 3. 

5. Application to a few AGN masing disks 

In the case where position-velocity diagrams are avail- 
able in the literature, we analyzed a few AGN masing 
disks to try to put constraints on masses. This is not 
a new problem regarding the mass of black hole. The 
mass of disks is, however, much a subject of debate, 
since there is almost no systematic or generic study like 
the one pr e sented here, except iso l ated attempts (see e.g . 
Hurell2002l; lLodato fc Bertinl 120031: iHerrnstein et all [20051 
Kondratko et al.ll2008l: iMcCallum et al.ll2009ll . Table 131 lists 
the systems considered here. Position-velocity data points 
were obtained by digitalizing graphs when published, with- 
out special treatment (in particular, we fully trust in the 
analysis by the authors to furnish deprojected position- 
velocity data in the reference frame of the black hole). 
The table also contains masses previously proposed for the 



black hole and disk, mostly through models of spherical 
distributions or models. In each case listed in Table [31 we 
have tried to fit the data by the linear formula and by the 
non linear formula as well. Position-dynamical mass dia- 
grams, together with the linear fit and averaged value, are 
displayed in Fig. @] The results are given in Table SJ In 
particular, from the linear regression, we estimated the de- 
gree of dispersion by computing the standard deviation a 
(column 8). From the simple analysis performed in Sect. 
HJ we immediately expect no precise values for NGC4945 
and Circinus. For NGC4258 and NGC1068, the dispersion 
is much weaker with a « 0.2 — 0.3. For these two objects, 
the linear approach should give the mass of the most mas- 
sive component quite correctly (i.e. see Tab. [5]). Finally, for 
IC1481, UGC3789 and NGC3393, dispersion is less than 
~ 10%, so both mass components should be obtained with 
"accuracy" . 

It was not possible to derive a triplet (Mbh, Mj, s) for 
all these systems through the non linear approach (negative 
mass or non convergence of the fitting procedure) , even for 
UGC3789 and NGC3393 whose data, like for IC1481, are 
spread over a large radial domain and little dispersed. In 
two cases (NGC4945 and Circinus), the method is inappro- 
priate: no correlation really exists in the position-dynamical 
mass and inferred parameters are non physical (in particu- 
lar, Md < for NGC4945). We are aware that some disks 
have been predicted to be geometrically thick, so the model 
of flat disk used here should naturally fail. In contrast, re- 
sults obtained through the linear approach are satisfactory 
for UGC3786, NGC3393, and NGC1068 where disk mass 
is similar to, or higher than the mass of the central black 
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Fig. 4. Position-dynamical mass diagram for each AGN listed in Table [3] ( left) and deviation in radius with the linear 
fitting (right). Observational data are circles, and the red dashed line is the averaged value of the {fli} considered. The 
linear fit through Eq. [TH]is the bold line (see Table H]). 



hole. For these objects, the massive disk should be prone 
to instabilities. For NGC4258, the slope is very weak and a 
mass ratio of about 0.05 is expected. 



6. Concluding remarks 

In this paper, we have reported a simple method to estimate 
the mass of a disk surrounding a central black hole and as- 
suming a flat, pressure-less disk at centrifugal balance. This 
model obviously applies to other kinds of systems contain- 
ing a disk and (possibly) a central star. Instead of the con- 
ventional position- velocity diagram, the position-dynamical 
mass diagram makes the measurement of possible devia- 
tions to Kepler's law easier. Using recent calculations for 
the gravitational potential of "truncated self-similar" disks, 
we have shown that i) the rotation law is generally not a 
power law of the radius, and ii) such deviations are, at zero 
order at least, directly related to the cumulative disk mass. 



It is clear that the disk model is very simple and can 
be improved in several ways. At the same time, it is dif- 
ficult to make the method robust and universal, since 
position-velocity data generally show a certain dispersion 
that can have several origins like thickness effects, instabil- 
ities (warps), etc. Through a simple analysis of uncertain- 
ties, we have shown that the linear approach gives quite 
correctly the most massive component (typically with a few 
tens of percent) if data dispersion is large (standard devi- 
ation larger than 0.1). For weak dispersion, both the disk 
mass and the black hole mass should be accessible through 
non linear data fitting, as well as the surface density power 
index. Moreover, the systemic velocity plays a major role, 
and the asymmetry often observed between the redshiftcd 
part and the blueshifted part of the rotation curve must 
be accounted for. This method must therefore be seen as 
a first step in the analysis of masses in star/disk systems 
based on gravity. 



8 



J.-M. Hure et al.: AGN disks and black holes on the weighting scales 



Acknowledgements. We thank S. Collin and F. Herpin for useful com- 
ments. The referees are acknowledged for their comments and sugges- 
tions, in particular about the problem of uncertainties. 



References 

Abramowicz, M. A., Czerny, B., Lasota, J. P., & Szuszkiewicz, E. 

1988, ApJ, 332, 646 
Binney, J. & Tremaine, S. 1987, Galactic dynamics (Princeton, NJ, 

Princeton University Press, 1987, 747 p.) 
Braatz, J. A., Reid, M. J., Greenhill, L. J., et al. 2009, in Astronomical 

Society of the Pacific Conference Series, Vol. 402, Astronomical 

Society of the Pacific Conference Series, ed. Y. Hagiwara, 

E. Fomalont, M. Tsuboi, & M. Yasuhiro, 274-+ 
Collin, S. & Zahn, J. 2008, A&A, 477, 419 
Collin-Souffrin, S. & Dumont, A. M. 1990, A&A, 229, 292 
Combes, F. 2001, in Advanced Lectures on the Starburst-AGN, ed. 

I. Aretxaga, D. Kunth, & Ft. Miijica, 223-+ 
Conway, J. T. 2000, MNRAS, 316, 540 
Dubrulle, B. 1992, A&A, 266, 592 
Evans, N. W. & Collett, J. L. 1993, MNRAS, 264, 353 
Evans, N. W. & Read, J. C. A. 1998, MNRAS, 300, 83 
Greenhill, L. J. & Gwinn, C. R. 1997, Ap&SS, 248, 261 
Greenhill, L. J., Moran, J. M., & Herrnstein, J. R. 1997, ApJ, 481, 

L23+ 

Herrnstein, J. R., Moran, J. M., Greenhill, L. J., & Trotter, A. S. 

2005, ApJ, 629, 719 
Heyvaerts, J. F. & Priest, E. R. 1989, A&A, 216, 230 
Hure, J. 2002, A&A, 395, L21 
Hure, J.-M. 1998, A&A, 337, 625 

Hure, J.-M., Hersant, F., Carreau, C, & Busset, J.-P. 2008, A&A, 
490, 477 

Kondratko, P. T., Greenhill, L. J., & Moran, J. M. 2008, ApJ, 678, 
87 

Lodato, G. & Bertin, G. 2003, A&A, 398, 517 

Mamyoda, K., Nakai, N., Yamauchi, A., Diamond, P., & Hure, J. 

2009, PASJ, 61, 1143 
McCallum, J. N., Ellingsen, S. P., Lovell, J. E. J., Phillips, C. J., & 

Reynolds, J. E. 2009, MNRAS, 392, 1339 
Mestel, L. 1963, MNRAS, 126, 553 

Miyoshi, M., Moran, J., Herrnstein, J., et al. 1995, Nature, 373, 127 
Papaloizou, J. C. B. & Lin, D. N. C. 1995, ARA&A, 33, 505 
Pringle, J. E. 1981, ARA&A, 19, 137 

Reid, M. J., Braatz, J. A., Condon, J. J., et al. 2009, ApJ, 695, 287 
Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 
Shlosman, I. & Begelman, M. C. 1987, Nature, 329, 810 
Shore, S. N. & White, R. L. 1982, ApJ, 256, 390 
Uscanga, L., Canto, J., & Raga, A. C. 2007, ApJ, 663, 857 
Yamauchi, A., Nakai, N., Sato, N., & Diamond, P. 2004, PASJ, 56, 
605 



Appendix A: Note on the power index of the 
rotation curve 

Let 7 be the local index of the rotation velocity, namely 

(A.l) 



where 



v = \/GM i? 7 , 

where Mq is a reference mass. In the present model, it is 
not a constant, but a function of the radius. We easily find 
from Eq. [5] 

-77-^ = 1 + 27, (A.2) 
amzu 

with 7 — — h for a Keplerian motion associated with a 
central point mass Mq . From Eq. |H1 we also deduce 



d In ji M<j 



P 2 (s)(s + 2)w 



s+2 



d In w fi 
and then 7 can be written in the form: 

1 



7 



<>7kep. , 



(A.3) 



(A.4) 



^7ko P . = ^7 F 2(s)(s + 2)w 



s+2 



(A.5) 



is the deviation to the Keplerian index. For s = —2, we have 
<^7kcp. = 0: the rotation curve resembles a Keplerian curve 
due to a point mass (see above). This result was already 
known. The mean deviation is 



(<57kep.) = 



1 



1 



' f 1 ^P2(s)(s + 2)m s+2 dTZJ, (A. 
Ja 



0) 



where q = M^/M-qu is the disk-to-black hole mass ratio. 
The above expression can be integrated exactly for some 
values of the s-parameter. For instance, with s = — 1, we 
find (c>7kcp.} ~ f- It means that a deviation to Kepler's 
law as low as 0.02 on the velocity index (i.e. 7 = —0.48) 
could imply a mass ratio q = 10%. This value is, in mag- 
nitude, consistent (within a factor ~ 2) with the monopole 
approximation that predicts a mass ratio ~ 0.4. For sys- 
tems containing a disk significantly less massive than the 
central object (i.e. q < 1) and A <C 1, we find 



(^7ko P .) 



2) 



1 + q 2{s + 3) 



P2(S). 



(A.7) 



Figure IA.1I displays the mean velocity index (7) versus 
q for a few values of the s-parameter. This plot enables ei- 
ther to predict the maximum departure to Kepler's law for 
a given mass or to bound the disk mass given a mean veloc- 
ity index. Practically, if deviations remain of small ampli- 
tude for low-mass disks, these may not be exploitable (data 
dispersion, thickness effects). 




0.01 """ 0.1 1 
mass ratio q 

Fig. A.l. Relation between the mean index of the super- 
Keplerian rotation law and the disk mass to black hole mass 
ratio for some values of the s-parameter. To a given mean 
index (7), it corresponds a minimum mass ratio q (red). To 
a given mass ratio, it corresponds a maximal mean index 
( blue). 



